% 创建时间：2025-06-21 12:51:24
% 作者：兪馳
clc;
clear;

A=[8,1,6;
    3,5,7;
    4,9,2];
Q=[16,4,1;
    9,3,1;
    4,2,1]

x=lyapsym(sym(A),-inv(A'),Q*inv(A')),
norm(A*A*A'-x+Q)